Critical dynamics of an isothermal non-ideal fluid 
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A pure fluid at its critical point shows dramatic slowing down in its dynamics, owing to a di- 
vergence of the order-parameter susceptibility and the coefficient of heat transport. Under isother- 
mal conditions, however, sound waves provide the only possible relaxation mechanism for order- 
parameter fluctuations. Here, we study the critical dynamics of an isothermal, compressible non- 
ideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydro- 
dynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics 
of an isothermal fluid effectively reduces to "model A", characterized by overdamped sound waves 
and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. 
Possible applications of the model to monolayer films are discussed. 

I. INTRODUCTION 

The presence of long-range static correlations can induce drastic slowing down of the dynamics of a fluid at the 
critical point [l|, Q • This can be understood by noting that the order-parameter relaxation rate is typically given 
as a ratio between a transport coefficient and a susceptibility: the susceptibility diverges critical point, while the 
transport coefficient remains roughly constant (or has a much weaker divergence). The characteristic dependence of 
the relaxation rate T on the wavenumbcr, r oc k z , or, by virtue of the dynamic scaling assumption, cquivalcntly on 
the correlation length, T cx £ -z , defines the dynamic critical exponent of the order parameter, z. In addition to the 
order-parameter relaxation time, other transport coefficients of a fluid, such as the viscosity, are often divergent as 
well and entail their own critical exponents. Similarly to statics, many properties of dynamic critical fluctuations, 
such as dynamic critical exponents or amplitude ratios, are universal and thus not specific to a particular substance. 
For instance, an energy-conserving pure (i.e., single-component) fluid at the liquid- vapor critical point and a binary 
fluid at the demixing point both belong to the same dynamic universality class of model H and thus share the same 
set of dynamic critical exponents @ . 

While the conventional pure and binary fluid have been extensively investigated both by theory and experiments 
(see, e.g., 04^] for reviews), the critical dynamics of an isothermal compressible, single-component fluid seems not to 
have received much attention so far. Two-dimensional isothermal fluids are often used as simple models for monolayer 
films that are confined to liquid interfaces 04^1 • Indeed, many of these films are also known to undergo liquid- vapor- 
like phase transitions |10l - [l2l ]. Recently, there has been growing interest in understanding the critical properties of 
these and related lipid lipid bilayer systems fl3l - |22j , especially, since they constitute the building blocks that form the 
membranes of biological cells [23[ . Since mono- and bilaycrs have critical temperatures of the order of the human body 
temperature, their critical dynamics might also have important consequences for physiology I2ll. |23I425{ . Isothermal 
non-ideal fluid models have also been used to study phase-separation [26h32| , capillary waves [33l |34| and supercooled 
liquids close to the glass transition (3^413 • All these works, however, did not address the critical dynamics of an 
isothermal fluid. 

In this work, using a mode-coupling approach, it is demonstrated that the isothermal condition leads to decisively 
different dynamic critical properties than in the standard model H universality class. A scaling analysis of the leading 
self-energy contributions emerging from the nonlinear Langcvin equations shows that the order-parameter dynamics 
effectively reduces to a time-dependent Ginzburg-Landau model for a non-conserved order parameter, known as 
model A Q. The upper critical dimension, d c = 4, is the same in statics and dynamics. The bulk viscosity diverges 
at the critical point with an exponent larger than in the mean-field limit, leading to overdamped sound modes at 
criticality. The shear viscosity, in contrast, remains finite, except for a possible weak divergence in two dimensions. 
The theoretical analysis is complemented by lattice Boltzmann simulations of the fluctuating hydrodynamics equations 
in two dimensions. We find that, even in this low dimensionality, the values of the critical exponents for the order 
parameter and bulk viscosity agree well with the analytical predictions that are obtained close to d c based on pure 
modcl-A behavior. 

In order to appreciate the difference of isothermal critical dynamics from that of a non-isothermal fluid, it is useful 
to recapitulate the basic results of the model H universality class. The original, incompressible model H consists of an 
advection-diffusion equation for the order parameter <f>, which is coupled to a transverse velocity field u [2Tltjl. l42ll50ll5lj |. 

d t (f>=~V-( ( f > u)+X\7 2 — .F+7T, (1) 
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z — a/v 
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TABLE I: Comparison of characteristic properties and critical indices of an isothermal compressible fluid and an energy- 
conserving pure fluid (or a binary fluid at the demixing point) described by model H. The cited numerical values are rounded, 
see the original works for more detailed predictions. Remarks: 1 "non-classical" refers to deviations from predictions of van- 
Hove theory, which assumes constant kinetic coefficients (see text). 2 Theoretical predictions (averaged) based on model A 
[iol. l4lj|. 3 Prediction of the scaling theory. Present simulations, however, only reveal a finite critical contribution to the shear 
viscosity. See text for further discussion. Extrapolation of the e-expansion results 2, S3] to 2D. 5 see [4314471 ]. 6 see [4£|. Static 
critical exponents have Ising values 0, Sj| and are identical for the isothermal and conventional (model H) fluid, v is the 
correlation length exponent, 7 = (2 — r\)v the susceptibility exponent, r\ the anomalous dimension exponent and a the specific 
heat exponent. 



Pod t u = - J^j +C s V 2 u + 7f ± , (2) 

where the label _L indicates that the transverse projection should be taken, i.e., the projection orthogonal to the 
wavevector in Fourier space (cf. sec. [II]). In the above equations, po is the mass density of the fluid, J 7 is a Ginzburg- 
Landau free energy functional, A is a bare kinetic coefficient (thermal conductivity), £ s is the shear viscosity and it and 
7r are appropriate noise sources. Importantly, in the case of a pure fluid at the liquid- vapor critical point, the relevant 
dynamical order parameter <f> is the entropy density, whereas for a binary fluid at the demixing point, cj> represents 
the concentration. The longitudinal part of the momentum density is neglected in the original model H, as the fluid 
is assumed to be incompressible, V • u = 0. This is a valid approximation at criticality, since thermal conduction (or, 
correspondingly, concentration diffusion in a binary fluid) proceeds on a much longer timescale than the propagation 
of sound waves. The latter process is therefore irrelevant to the dynamics of the order parameter [Hoi. [52j. In model 
H, the order-parameter relaxation rate is given in Fourier space by 

with x(k) being the susceptibility (isothermal compressibility) This relation can be derived from eqs. ([1]) and © 

by linearizing and neglecting the advection term. Since x(k) oc k~ 2+n , the "classical" (van Hove [HH) theory predicts a 
dynamic critical exponent of z = 4 — 77 for model H. The classical result, however, turns out to be violated in a real fluid, 
since the kinetic or transport coefficients are, due to the presence of reversible mode- couplings, affected by the critical 
order-parameter fluctuations as well j2l. |54H56| . Specifically, in model H, the kinetic coefficient A is renormalized by the 
nonlinear advective coupling between cb and u, changing the dynamical exponent to z = 4 — 77 — z\ = d + y, where z\ 
and y are the exponents characterising the divergence of A and the shear viscosity ( s , respectively (4214441 [5lT l54l |57| . 
This value of z is also close to the predictions of early mode-coupling approaches, z ~ d [5(3, H3- The effect of critical 
fluctuations on the shear viscosity, £ s oc £, y , is comparatively weak, leading only to a small exponent y ~ 0.07 in 3D 
HH S3 HH H3 1 as confirmed by experiments [i^, S(| . 
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Critical fluctuations in a pure fluid have an effect on sound waves as well, which, however, is one-sided since the latter 
are decoupled from the order-parameter dynamics. Acoustic effects can be studied with an extended, c omp ressible 
version of model H that includes the full set of equations for the mass, momentum and energy density j^. l59l - l67| (see 
also (68j and references therein). In an energy-conserving pure fluid, sound waves propagate with the adiabatic speed 
of sound |6§{n|, 



dp 



r 2 

S ' ad dp 



= ^L-, (4) 
s c v px 



where p is the pressure and c p , cy are the specific heats at constant pressure and volume. These are related by 
c p = cy + Tx(3y/p, where /3y — {dp/dT)\ p is the thermal pressure coefficient (slope of the p — T-curve). Since /3y 
is not critical, we have c p ~ \ and hence the critical behavior of the speed of sound (at zero frequency) is given by 



oc r a/v ■ (5) 



For comparison, in a binary fluid, the critical sound speed is governed by the constant-pressure specific heat, whose 
divergence is - due to a larger background contribution - much less pronounced than for the pure fluid. The critical 
sound damping is given by the bulk viscosity, which can be determined from a Green-Kubo relation involving the 
nonlinear pressure fluctuations [48l [59l - l6l| . The extended model H predicts a strongly diverging bulk viscosity, £j oc £ x 
(at zero frequency), with x = z — a/v being ~ 2.9 in 3D. 

While theoretical and experimental investigations of critical dynamics in pure or binary fluids have a long history, 
simulations are scarce and have been performed only quite recently [7214771 ] . However, most of these studies are 
in 3D, and values of the dynamic critical exponents for model H in 2D seem so far not to have been obtained 
either by experiment or simulations (cf. [2ll [78l - [83j ). The values given in Tab. Ufor model H therefore correspond 
to extrapolations of theoretical e-expansion results @, |42[ to 2D (throughout this paper e = 4 — d). The shear- 
viscosity exponent y in Tab. U has been computed using the e-cxpansion result for z\ in conjunction with the relation 
y = 4 — d — r\ — z\, giving the lower bound, as well as the direct e-expansion result for y, giving the upper bound. 
Since the 0(e 2 )-term in the e-expansion of z\ is quite small relative to the leading term, one might suspect that the 
extrapolated value will not be grossly unrealistic. Of course, one has to keep in mind that, even in the non-critical 
case, transport coefficients in a two-dimensional fluid acquire logarithmic divergences in the long-time or -wavelength 
limit |84l - l87| . which might interfere with the critical divergences. 

From relation (jj) we see that a description in terms of the isothermal speed of sound, 



dp 
dp 



— , (6) 
PX 



would become applicable if the specific heat ratio c p /cy would be close to one. Dynamically, adiabatic conditions are 
achieved if the thermal relaxation rate [eq. ([3])] is much smaller than the characteristic frequency of a sound wave, 
i.e., 

r < c s k . (7) 

Since c s , a dfc oc j^ 01 /"- 1 at k ~ above relation is clearly fulfilled in a ordinary fluid close to the critical point. 
Relation (J7J) provides a posteriori also a justification for neglecting the "faster" acoustic processes in usual model H 
calculations. Far from criticality, violations of condition ([7]) can occur at finite wavenumbers in the hydrodynamic 
regime. Isothermal conditions can thus be only achieved by coupling to fluid to some kind of heat bath, so that 
temperature fluctuations are removed at a sufficiently fast rate. On the other hand, the friction between fluid and 
substrate should be sufficiently small, in order not to break momentum conservation and violate the characteristic 
sound mode behavior of the compressible fluid in the relevant wavenumber regime [88| . After these introductory 
remarks on the critical dynamics of ordinary non-ideal fluids, we now turn to the analysis of the isothermal non-ideal 
fluid. 



II. THEORY 



An isothermal compressible fluid is governed by a continuity equation for the density p and a conservation equation 
for the momentum j = pu [y, |26|, |33|, |69( , 



(8) 
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da 



-pv^- + c,v 2 - + (a 



CJl-2/d])VV-- + V-R. 



.SjF 

6p ' p 

Here, £ s and <^b are the bare shear and bulk viscosity and R is a random stress tensor with correlations [69l |8S 
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(R a p(T,t)R yS (r',t'))=2k B T 



5{r - r')S(t - t') 



(9) 



(10) 



imparting Gaussian thermal noise on the fluid. The static probability distribution of the density fluctuations are 
governed by the Ginzburg-Landau free energy functional 



T 



I 



dV 



^(V0) 2 



/oW>) 



where n is a constant and 



i> = {p- po)/po 



(11) 



(12) 



is the order parameter, po is a constant background density. The Landau free energy density fo is taken to be a 
quartic polynomial in (j), 



/o(0) - \r<f + . 

The "streaming term" involving J- in eq. ([9]) can be written as divergence of a pressure tensor P [33l . |90j-l9 

gjr K 

pV— = V • P = Vp - — pVV 2 p 
Sp p z 



where 



P. 



a0 



[PO 



and the scalar pressure po is given by 



Pa = pdpfa - fo 



r<f> + -r<f> + u<jf + 



(13) 



(14) 



(15) 



(16) 



An essential complication in the analysis of the compressible Navier-Stokes equations is the presence of the non- 
linearity j/p in the viscous stress (36l . |39i . [95| . Here, we treat this term perturbatively by expanding 1/p around the 
background density po, i.e. 1/p = 1/po — (l/pfySp + . . . [§6|, We find that this term is irrelevant for the critical 
behavior of the fluid above two dimensions. Furthermore, in the above Navier-Stokes equations, the convection term 
V(puu) has been neglected, since the effect of this term is well known and not specifically related to critical phenom- 
ena [84148 7j . Above two dimensions, the convective nonlinearity gives a finite contribution to the renormalization of 
the shear and bulk viscosity, while logarithmic divergences appear in 2D in the limit of large time- or length-scales. In 
principle, these effects require careful treatment, also regarding a possible interplay with critical fluctuations. This, 
however, is out of the scope of the present work, and in the subsequent analysis, it is therefore assumed that these 
effects are either already included in the bare viscosities or can be considered separately from critical fluctuations. 

To proceed, eqs. ([SJ and (O are written in Fourier space, 



ujdp = k • j , 



(17) 



- iwj = -ikc 2 s (k)Sp - ikp„i - iN - v s k 2 [j + Y] - {y h +v s [l- 2/rf])kk • [j + Y] + ik • R + ik ■ h , (18) 

where v s = ( s /po and v^, = (b/po are the kinematic shear and bulk viscosities. The generalized isothermal speed of 
sound [Hi 

c 2 s {k) = c 2 s + K Po k 2 =r/p + np k 2 , (19) 

contains the linear part of the thermodynamic pressure, while p n \ and N arc the Fourier-transforms of the remaining 
nonlinear parts: 

Pnl (k, u) = \rl <t>{k - q)m +uf_ Hk-q- ?)m<K?) + |« / - ? - - f'MffltfM?') , (20) 

J q J q,q' J q,q' ,q" 



5 




-2 -1 1 2 -4 -2 2 4 



(a) (b) oj/(2nT) 

FIG. 1: Typical shape of the dynamic structure factor C(k, ui) in an isothermal fluid in the weak-damping (a) and strong- 
damping (b) regime. 



N(k,u) = K / <K% (21) 

Jq 

Here, the shorthand notation k = (k, u), q= (q, a), etc. and J. = J is introduced. The quantity 

Y(k,w) = lm<f>{k-q) (22) 

Jq 

represents the leading correction term of the expansion of j/p around j/po- 

The nonlinear Navier-Stokes equations (fTT)) and (fT5)l can be split into longitudinal and transverse parts with respect 
to the wavevector k. The corresponding longitudinal (T) and transverse ('t') projections of a vectorial quantity 
v = v/k + v t are defined as vi = k • v and v t = 7kV, where k = k/fc and 7k = (I — kk). Analogously, for a tensorial 
quantity R, we have k • R = i?;k + Rt, with Ri = k • R • k, and Ht = k • R • 7k- We thus arrive at 

uiSp = kji , (23) 
cvji = kc 2 s (k)Sp + k Pnl + N[- mk 2 (ji +Y t )-kRi- kh t , (24) 

wjt = N t - mk 2 (it + Y t ) - kRt - kh t , (25) 

where i>i = v j, + i^(2 — 2/d) denotes the longitudinal and v t = v s the transverse viscosity. The longitudinal and 
transverse parts of the random stress tensor are correlated as (|i?;(k,a;)| 2 ) = 2 po^ik bT and (|i? tjCt (k,a;)| 2 ) = 2p v t kBT. 
Combining eqs. (|23|) and ([2~4")l . the longitudinal current j; can be eliminated completely, leaving only eq. (f2"5")) for the 
transverse current and a single, nonlinear sound-wave equation for the order parameter (setting henceforth po = 1): 

- co 2 (j) + k 2 c 2 {k)<j) - muk 2 c/) = -k 2 Pn i - kNi + iiyik 3 Yi + k 2 R l + k 2 h t . (26) 

Obviously, the term Y, which can be written as 

H~k - q) , (27) 

provides a bi-directional coupling between the order parameter and the transverse current. Additionally, the transverse 
current is affected by the order-parameter fluctuations through the term N t . Note that, due to the way the free energy 
functional enters the Navier-Stokes equations, there appear to be more nonlinearitics in eq. (|26[) than the corresponding 
static critical theory or in nonlinear sound-wave equations studied in isotropic elastic phase transitions (98l . |99| . In 
the latter, the restoring force is given by a term ~ 5T /5p rather than by a pressure gradient. 



■h(q) 



A. Linear hydrodynamics 

At this point it is useful to collect some basic results of the linear model. Neglecting the nonlinear terms in eq. (|2GI) . 
the bare response- and correlation function (labeled by an index 0) for the order parameter are given by 

k 2 
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[w 2 - £; 2 c 2 (k)] + (uv l k 2 ) 2 
The response and correlation functions are related by a fluctuation-dissipation theorem, 

Oh. 

Co(k,w) = — — ImG (k,w), (30) 
which, in the zero-frequency limit, becomes the fluctuation-response relation 

C (k) = k B TG o (k,u = 0), (31) 

where C(k) is the static structure factor, 

<7 (k)= / ^C (k )W ) = ^. (32) 
7 2tt c 2 (k) 

In the linear case, cq. (|26|) represents a damped harmonic oscillator driven by random noise [49| . The dispersion 
relation of the associated sound waves is given by 

u = ± [c^(fc)fc 2 - vfk A /A] 1/2 - rf/2 • (33) 
In the weakly damped case, where c 2 (fc) > vfk , sound waves are oscillating with frequencies uj a = 

1 /2 

±/c [c 2 (fc) — i/j 2 fc 2 /4] w c s (k)k and are exponentially damped with a rate of vik 2 /2. In the opposite, overdamped 
case, the solution (|33|l becomes purely imaginary and sound waves decay, in the limit of long times, with a rate of 

r(k) = . (34) 

VlR 

At short times, there is another decay regime with a rate ~ fik 2 , which, however, is negligible if c 2 3> ifk . The 
response and correlation functions in the overdamped, long-time limit can obtained by neglecting the "inertial" term 
w 2 in eqs. ([25]) and flM), yielding 

G (k,w) = — 1 — — jz-r , (35) 

C (k,w) = — 77T - 2~TT^7TT ' ( 36 ) 
viR\k) u z + r-(k) 

In the time-domain, this corresponds to a pure exponential decay: 

G (k,i) = — !— caq>[-r(k)t]0(t), (37) 
viR\k) 

C (k,i)=C(k)exp[-r(k)|i|]. (38) 
Similarly, for the transverse current, we obtain from cq. (|25[) the bare response and correlation functions 

G tfi (k,u)= * (39) 

U! + lft k 

2v t k B Tk 2 

Ct,o(k,a;)= 40 
w 2 + [ytk 2 -) 1 

where (7t a (k, uS)jt t p{k! ,u)')) = Ct.o(k, uj)(2ir) d+1 5(k+k')5(uj+uj')8 a i3. The static correlations of the transverse current 
are independent of the wavenumber, 

C t (k) = k B T. (41) 
The linear hydrodynamics expressions (f28|) and (|29|) can be cast into standard dynamical scaling forms [U, |49|, llOC| , 

G (k,c) = e^Gi^^n, c (k,u) = e~ v+z c{k£,M z ) , m 

with Q and C being scaling functions, r/ = and a dynamic scaling exponent of z = 2. This value of z can also be 
inferred from the damping rate in the overdamped case, eq. (|34[) . Analogously, for the transverse current, eqs. (|39[) 
and (|40p . we have 

G t ,o(k, u ) = f- 1 g i (fe^f) ! C t ,oCk,w)=e*'C t (fce,a;e**) J (43) 
with a dynamic exponent of z t = 2. 
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FIG. 2: Fundamental vertices of the model arising from eqs. (|20[) . (|21|l and (J22J . Here, k is an external wavevector, while q and 
<r denote an internal wavevector and frequency. A filled circle represents a coupling constant and an integration over internal 
wavevectors and frequencies respecting space- and time-translational invariance. 
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FIG. 3: Leading frequency-dependent diagrams contributing to the self-energy of the order parameter, E(k, u). A solid (wavy) 
line with an arrow represents Go (Gt,o) an d a solid line with a open circle Go (Ct,o)- Dashed lines indicate amputated external 
"legs" for clarity. 



B. Critical order-parameter dynamics 



The critical dynamics of the o rder parameter, governed by eq. is discussed here within a mode-coupling 

approach 0, [5(| [HI, [H, llOll 1 1021 ] . To this end we shall cons truct a p e rtur bative solution of the nonlinear order- 
parameter equation using the response function formalism (55l . llOOl . Il03l Il04| and identify the leading contributions 
via a scaling analysis. With the help of the bare response function Go, eq. ([26]) can be rearranged as 

(44) 



(j) = fa - G p„i - G Ni/k + mkG Yi + G hi , 



with 0o = Go-R; being the zeroth-ordcr solution. Due to the coupling between the order parameter and the transverse 
current, we also need to consider eq. (|25j) . which can be written as 

(45) 



jt = -jt,o + G ti0 N t /fc - mkY t - Gt,oht 



with j t , = Gt )0 Rt. 

The nonlinearities on the right hand side of eqs. (|4"4"]) and (|43|) can be translated into a diagrammatic representation 
as given by Fig. [2] There, a solid (wavy) line with an arrow represents Go (Gt,o), a thick solid line represents the order 
parameter (f>, a filled circle stands for a coupling constant and an integration over internal wavevectors and frequencies 
respecting space- and time-translational invariance. The vertices involving the couplings r, k, V\ and v t have no 
counterparts in the static theory or in standard Ginzburg-Landau models Q and arc specific to the compressible 
fluid. They arise already from the Gaussian part of the free energy functional and are known to be important in the 
case of a supercooled liquid [H, [39| . 

A solution of eq. (|44|) for 4> can be iteratively constructed, following standard rules [Hoi. IB"!! . llOOl . 1 1031 ] . and leads to 
a Dyson relation for the full response function G = 8(<f))/8hf. 



G(k,w) = 



Gq 1 (k, u) - S(k, u) -lo 2 + k 2 c? s (k) - mujk 2 - fc 2 E(k, u) 



(46) 



Here, £ is a self-energy, which encapsulates the effect of the nonlinear interactions between the order-parameter 
fluctuations. These can be understood to lead to a renormalization of the transport coefficients of the fluid of the 
form 



C yk)= C 2(k)-Re£(k,0), 



(47) 
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v m (k) = vi +n(k,w 0) 



(48) 



where II(k, w) = 9ImE(k, uS)/du). We shall focus here only on the small-frequency limit and neglect any frequency- 
dependence of c s r and vir. Thus, we assume that the response function keeps the same principle form as in linear 
hydrodynamics, but with appropriately renormalized transport coefficients. Also, a possible renormalization of the 
background density is neglected here. 

The nonlincarities in the model give rise to a large number of diagrams contributing to the static and dynamic 
(i.e., frequency- independent and -dependent) parts of the self-e nergy. T he dominant contribution can in principle be 
determined through a straightforward scaling analysis (see, e.g., [5J, ll02( ) , making use of the dynamic scaling forms of 
the response and correlation functions stated above. Regarding the static parts, we can just invoke the fact that the 
nonlinear Langevin equations of the model preserve the equilibrium probability distribution of the Ginzburg-Landau 
free-energy functional @, |89[. Thus, in the hydrodynamic limit, the renormalization of the isothermal speed of sound 
must be consistent with the static theory, implying that, c 2 = 1/px <x C™ 7 ^, where x 1S the isothermal com press ibility, 
£ the correlation length and 7 and v arc the usual static critical exponents [note 7 = (2 — if)v] 0, 0, Q |128| . The 
upper critical dimension for statics is 4. In contrast, in a conventional fluid (model H), sound waves propagate with 



the adiabatic speed of sound, which vanishes much more weakly at the critical point, c s oc £ 



-a/ ' v 



with a being the 



specific heat exponent [l|, H, |48|, [59|, |64J . 

Turning to the dynamics, the leading irreducible, frequency-dependent diagrams contributing to E are shown in 
Fig. El For a given diagram E(j), we have ~ £(i)£ z as far as the scaling behavior is concerned. We do not consider 
here an expansion in the number of loops, but rather focus only on the leading diagrams arising from each vertex. 
The analytic expressions of the individual diagrams are given by 



E„(k, w) 
E„,(k,w) = 54u 2 



18u z / G (k - q - q')C (q)C (q') 



G (k - q - q" - q")C (q)C (q')C (q") 



9, <?',<?" 



E r (k,w) 



+ 8lu 2 / G Q (k - q)C (q) 

Jq 

-r 2 f CoCk - q)Co(q) , 

Jq 



Co(q) 



E*(k,aO = n 2 / G (fc - q)Co(q) [(k ■ q) 2 + 0(k 4 )] , 

Jq 

E„j(k,w) - -vf f G (~k - q)C (q) [w<r(k ■ q) 2 + 0(uj 2 ) 



(49) 
(50) 

(51) 
(52) 
(53) 
(54) 
(55) 
(56) 



In eqs. (|53[) to (|56p . only the principle form of the kernels is indicated, which is sufficient to derive scaling properties. 
Also, expressions for the remaining one-loop diagrams of Fig. [3] that involve two different couplings or a transverse 
current response/correlation function are not stated explicitly but can be easily obtained. In fact, it will not be 
necessary to compute them explicitly, since all vertices involving r, K, v\ or v t scale in the same way, up to differences 
of 0(77). To see this, note that jt scales like </>£ -1 , as can be inferred from the form of the corresponding correlation 
functions [eqs. (J42J , (|43|) ]. Some of the additional diagrams at two-loop order are briefly discussed in appendix [A] 
they will give rise to subdominant contributions and can thus be safely neglected. 

First of all, we consider the mean-field approximation, where the values of the dynamic scaling exponents are given 
by z = z t = 2. Taking into account the strong temperature-dependence of the Landau parameter, r ~ 1/x ~ £~ 2 , we 
find, in the limit k — > cj — > £ -z : 



Ej(k, u) = v x v t I G t ,o(q)C (k - q) ujq (k • Tq ■ kj + 0(u)k) 
Go(q)C tfi (~k - q) [k-q + 0(fc 2 )] + ... . 



n u oc £ 



8-2d 



Tl v i oc £ 



10-3d 



n, 



,3 K C 



2-d 



(57) 



The identical scaling of all one-loop diagrams of Fig. [3] is a consequence of the identical scaling behavior of the two- 
point vertices in the model. We can also obtain scaling predictions beyond the mean-field case by making use of our 
knowledge of the proper critical behavior of the static couplings, k, r and u. These are renormalized by fluctuations 



9 



as r R ~ £,' 2+, \ k r ~ f, u R ~ £<*- 4 +2'? @,[lOO|, implying that 

II,, a £- a +" , II,,',,,, (58) 

The remaining contributions due to H, all scale oc £ z - d ±°(v) ^ or £ z t-d±0(ri) ^ respectively, up to differences in the 
exponent of 0(rj) accounting for possible divergences of v\ or v t . Thus, all diagrams except II U are irrelevant for 
d > 2, provided that z and Zt are still close to 2, which will indeed be the case. Note that the contributions from the 
u<j) -vertex, which were found to diverge below a critical dimension of d = 10/3 in the mean-field limit, now remain 
finite at least down to three dimensions, owing to the renormalization of u. We also see that the j/p-nonlinearity in 
eq. ((9]), responsible for the coupling between longitudinal and transverse current, is not relevant for d > 2 and it is 
safe to approximate p by po, as far as asymptotic critical properties are concerned. 

From the dominance of E u , which arises from the 3 -vertex of p n \, we conclude that the upper critical dimension of 
the present isothermal non-ideal fluid model is d c = 4 both in statics and dynamics. The present analysis also reveals 
that the relevant nonlinearities responsible for the deviations from the classical (van Hove) predictions are different 
for the isothermal fluid and model H: in the latter case, the deviation is caused by the reversible advection term [see 
eq. ((TJ)], whereas in the isothermal fluid, it is caused by the dissipative 4 -nonlinearity of the Ginzburg-Landau free 
energy functional. Hence, in the isothermal fluid, all dynamic critical effects are induced purely by thermodynamic 
quantities. 

The scaling result for n„ of eq. |58|) is not sufficient to determine the precise value of z. This can be done via a 
renormalization group (RG) calculation slightly below four dimensions. To this end, the wavenumber integrations in 
n„ of eq. (|49]) are performed incrementally in a shell Aoe _s < q, q' < Aq, where Ao is a cutoff and s denotes the RG 
flow parameter. The contribution linear in s, which we shall write as A(y{) vi, has been calculated in [98|, with the 
essential result that A(yi — > oo) = 61n(4/3)r7 and A(vi) > A(oo) for any finite vi. These estimates have been obtained 
at 0(e 2 ) in an e-expansion. The RG equation for the longitudinal viscosity then reads 

d s vi(s) = A( Vl )vi{s), (59) 

from which one concludes that, for any positive bare vi(0), vi(s) will grow along t he R G flow and asymptotically scale 
as e A (°°) s . Thus, in the hydrodynamic limit, which is reached for e s ~ Ao£ P. Il05| . the renormalized longitudinal 
viscosity behaves as 

viR oc viC , (60) 

with a critical index 

sc = 61n(4/3)77~1.7»?. (61) 

In the critical regime (fc ^> ^ _1 ), we have accordingly, vir cx i>ik~ x ' . With this result, we can show that sound waves 
must be overdamped in the critical isothermal fluid: using the fact that c 2 sR oc £~ 2+,; , the linear hydrodynamical 
condition for strong damping, c S R{k) <C kviR , becomes £~ 1+,? / 2 /fc + const, x k~ n l 2 <C £ x - Thus, for wavenumbers of 
order ~ we have £ v £ 2x , which is always fulfilled in the asymptotic critical regime since 2x > 77. Of course, 
we could also have kept, for d close to d Cl only the dominant u(/) 3 -nonlinearity in eq. (1261) and thereby recover the type 
of sound-wave equation studied in the context of isotropic elastic phase transitions [98| . The associated RG analysis 
in [98| then leads to the same predictions as above. 

Concluding, in the critical regime, eq. (|26p reduces in the long-time limit to model A in the classification of 0, Il06l - 

Il08l |. that is, a time-dependent Ginzburg-Landau model for a non-conserved order parameter, 

1 

- \oj5p = — + K + hi , (62) 

n op 

where 1Z = Ri/vi is a noise source of variance ~ \/vi [129| . Since overdamped sound waves relax with a rate 

T = ^ , (63) 

ViR 

the dynamic critical index z, defined via r ~ £ z , of the fully nonlinear fluid model follows as 

z = 2 - T] + x . (64) 

In contrast, in the linear hydrodynamic (mean- field) case, we have z = 2 and x = 0. If we assume that the pure model- 
A behavior of the critical isothermal fluid persists also in low dimensions, we expect, in the interesting two-dimensional 
case, a value of 



z~ 2.08. ..2.17 (2D) 



(65) 
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based on recent theoretical calculations and Monte Carlo simulations of m odel A liol l4ll Il09l | . In the wider literature, 
varying estimates for z, ranging between 2.0 to 2.3, have been reported jllOl . Ill I] . Above value for z translates to 
x ~ 0.4 and agrees surprisingly well with the O(e 2 )-renormalizati on-g roup prediction of eq. (j6Tj) in 2D. For comparison, 
for a conventional fluid (model H), z ~ d and x = z — a/u [1, 0, l48j 

Returning to the scaling estimates of eq. (|58[) , a value of z > 2 would imply the weak divergence of various diagrams 
in 2D, which could provide corrections to the critical exponents. To address this issue, explicit calculations of the 
corresponding contributions will be required. Our simulations (sec sec. IIIip . however, yield a value of z w 2.2 ± 0.1, 
suggesting that possible corrections to the model-A behavior are small at least. 

C. Critical shear viscosity 

The shear viscosity is discussed in the following based on a Green-Kubo approach. We consider the x-component 
of the nonlinear NSE and choose the wavevector to lie along the y-dircction, i.e. k = (0, k) in 2D. Applying the 
approximation j/p ~ j/po + (i/ p^)8p to the viscous stress tensor, whose xy-component becomes 

Sxy = v s [d x (j y + 3 y 4>) + d y (j x + j x 4>)} (66) 

the equation for the transverse current can be written as 

d t j x = ~v s k 2 j x - ikP xy (k) + ikSf y in {k) + ikR xy , (67) 

where P xy (k.,t) = —n J q x q y c/)(q)0(k — q) is the Fourier-transform of the off-diagonal term of the thermodynamic 
pressure tensor [eq. ([To])], while ^"^"(k, t) = iu s k f q ^(q, £)</>(k — q, t) contains the nonlinear terms of S xy [eq. ([66]) ] 
involving the order parameter. 

The fluc tuation contribution to the shear viscosity, f s , C rit) can now be inferred by invoking a Green-Kubo relation 
, ItH Ill2l | . For the contribution from the thermodynamic pressure tensor we find 

where, as usual, the four-point correlation has been decoupled into products of two-point correlation functions. For 
the contribution due to the nonlinear part of the viscous stress tensor, S xy , one writes j x (q) — jii^Qx + jt,x(*l) and 
uses the fact that the correlation function of the longitudinal current fulfills Ca(q, t) = dfC(q,t)/q 2 . Also, ji and jt 
are independent to leading order. This gives analogously 

1 r°° 

<crit = TTTT^F J™ J dt (^ in (k, t)S^(-k, 0)) oc ^-d-^v + f t -v-d , (69) 
Vk B T fe^-i Jo 

where the external wavevector is taken at £ -1 and, for the evaluation of the part involving the transverse current, 
it has been assumed that z t < z. In the mean-field limit, all contributions scale oc £ 2 ~ d and thus are finite for 
d > 2, while a potential logarithmic divergence is indicated in 2D. If, in contrast, scaling exponents appropriate for 
the true critical point are taken (where z ~ 2.2 and kr oc ^ ?? ), the contribution to ^ SiC rit of eq. (|68[) attains a weak 
power-law divergence in 2D, characterized by a critical exponent y = z — 2, implying that z t = 2 — y < 2 in 2D. As 
a consequence, eq. (|69[) becomes now finite in all dimensions. In our simulations, however, where we evaluate the 
critical shear viscosity using eq. (|68[) . we do not observe a divergence in 2D. This might be related to limitations of 
the simulation model and is discussed further in sec. IIIIB 21 Note that the same calculation as in eq. (f68|) applies to 
model H as well. There, the shear viscosity diverges mildly in 3D, due to a value of z ~ 3.07 that is slightly larger than 
3 0,S3]. In 2D, the uncertainty in the theoretical value for z (see Tab.|I]) and the lack of experimental or numerical 
studies permits no definite conclusion on a possible critical divergence of the shear viscosity in an ordinary fluid. 

III. SIMULATIONS 

The theoretical predictions are now compared to full fluctua ting hydrodynamics simulations of a non-ideal fluid 
in 2D, using the Lattice Boltzmann (LB) model introduced in 1 1 1 3| . Resu lts on the static critical behavior of the 
model, which belongs to the 2D Ising universality class, can be found in [l!4| . The two-dimensional case is interesting 
for several reasons: first, the isothermal condition could probably be realized here most easily experimentally by 
coupling the fluid to a heat-absorbing substrate. Second, the scaling arguments of sec. HI1 suggest that the fluctuation 
contributions of various nonlinearities grow around a dimension of d = 2. As such effects are difficult to assess 
analytically, numerical simulations can provide useful insights and are complementary in this case. 
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A. Setup 

Parameters of our LB simulations are chosen as in Ill4j . a typical setup at the critical point being r = —4.8 x 10~ 5 , 
u = 2.8 x 10~ 2 and k = 9.6 x 10 -5 . This choice le ads t o a m ean-field interface width of ~ 2 l.u. and is expected to 
ensure reliable results on the fluctuation dynamics [ll4l Ill5| . The noise temperature is set to kgT = 10 -7 and the 
bare bulk and shear viscosities to vi, = v s = 0.04/3. Quantitatively similar results have been obtained also for other 
parameter combinations. Simulation boxes are of size L x L = 256 2 , except for Figs. [SJandlBfc, where L x L = 128 2 . 
All simulation results reported in the paper are obtained with a standard LB implementation, where the viscous stress 
consists of terms of the form vpd a U/3, i.e., the dynamic viscosities depend on p. In a few cases it has been checked, by 
using an implementation where the p-dependence of the dynamic viscosities is eliminated, that results are not affected 
by this LB specific peculiarity. Besides, due to requirements of numerical stability, relative density fluctuations 5p/po 
in our simulations are on average well below a few percent, thus an approximation p ~ po is warranted. Due to 
the multiplicative nature of the renormalization of the relaxation rate, the specific values of the viscosities are not 
important. They do, however, influence the extension of the overdamped acoustic regime and the crossover from 
mean-field to the expected model-A critical behavior. While a small longitudinal viscosity leads to a more rapid 
equilibration of the order parameter, it also shifts the onset of the overdamped regime to smaller wavenumbers. In 
consequence, larger simulation boxes would be required to reduce the residual speed of sound at the critical point 
sufficiently. Furthermore, a large value of c s can significantly affect the long-time decay of the order-parameter 
correlation function, which can be misinterpreted as caused by a larger viscosity. In order to diminish these and other 
undesired finite-size effects, the lowe st fc-modes are usually excluded from the analysis of our results. To eliminate 
effects of lattice anisotropy (cf. [ll6l |). all wavenumber-dependent quantities shown in the plots are computed as an 
average over the Cartesian axes of the Fourier plane. 

The theoretically predicted logarithmic divergence of the viscosities in 2D due to the convective nonlinearity [84T - [87| 
is difficult to observe and requires either very long simulation times or large simulation boxes (cf. [113). Indeed, since 
the effect is proportional to fcgTlogL [Hf|, it is expected to be below the threshold of statistical accuracy for the 
present setup. 



B. Results 



1. Order parameter 

Fig. 0H shows the dynamic structure factor C(k,t) at the critical point for different wavenumbers k as obtained 
from our simulations. The exponential decay of C(k,t) is clearly seen in the semi-logarithmic plot. According to the 
dynamic scaling hypothesis, C(k,t) = fc _2+?? C((fc^)~ 1 , k z t); hence, sufficiently close to the critical point, where (fc£) _1 
is small, the dynamic critical index z can be determined by plotting C(k,t)/C(k,0) versus the rescaled time k z t, 
testing different values of z until a complete data collapse is achieved. This is done in Fig. 0t>, from which we infer 
a value of z ~ 2.25 ± 0.1. For comparison, the insets demonstrate that, when rescaling the data with a significantly 
larger or smaller value of z, the data collapse remains incomplete. 

Alternatively to the rescaling procedure, the dynamic critical index can be more directly determined from the 
relaxation rate which can be obtained by fitting an exponential decay [eq. (|38[) ] to the dynamic structure factor. 
The assumption of an exponential relaxation is well satisfied in the overdamped regime, after neglecting the short- 
time, non-exponential part of C(k, t) caused by a finite residual speed of sound. In Fig.|3t, the so obtained relaxation 
rate is plotted against the wavenumber in a double-logarithmic representation. At small wavenumbers, the expected 
power-law behavior r cx k z is clearly seen. A numerical fit yields a value of the exponent of z ~ 2.2 ± 0.2, which 
agrees well with the value obtained above by rescaling the structure factor data. 

Within the error margins, the value of z extracted from the simulations thus clearly shows that the dynamic critical 
index is significantly increased over its mean-field (z = 2) or van Hove (z = 1.75) value. Moreover, it agrees reasonably 
well with the theoretical predictions based on pure model- A-type critical behavior [eq. (|6"S")) ]. Note that the obtained 
value of z in principle also depends on the range of wavenumbers considered. The deviation from a pure po wer-law at 
larger wavenumbers might be caused by the general wavenumber-dependence of the LB transport coefficients |ll8l . lll9j | 
or the influence of other nonlinearities in the model that overwhelm t he leadin g order critical divergences. For k > 1 
also the discrete nature of the lattice becomes noticeable (see, e.g., jl!3l Ill4( ). More precise values for z could be 
obtained by using larger simulation boxes, thereby extending the low-fc regime and decreasing the influence of the 
non-exponential decay of C(k,t) at short times. 

In Fig. [5J the behavior of the order-parameter relaxation rate (obtained from exponential fits) is investigated in 
greater detail for different wavenumbers k and reduced temperatures 9, where 9 = (r c — r)/r c . As seen in Fig. EH, the 
relaxation rate at fixed k first markedly decreases when approaching the critical point (8 = 0), but eventually levels 
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FIG. 4: (a,b) Dynamic structure factor at the critical point: (a) Raw data for the lowest wavenumbers k (increasing from top 
to bottom). The time is given in lattice units and the lines are drawn as a guide to the eye. (b) Test of the dynamic scaling 
form of C(k,t). The value of the dynamic critical index z is determined by varying z until all points fall onto a master curve. 
The inset shows that the collapse is incomplete for z significantly different from 2.25. (c) Relaxation rate F at the critical point 
obtained from exponential fits to the dynamic structure factor for different wavenumbers. The solid line is oc fc 2 ' 2 . 




off at a finite value caused by a nonzero speed of sound. Also, the expected temperature dependence T oc 9 ZV with 
zv ~ 2.2 is not observed, but instead a less steep decrease. Th ese e ffects are well known consequences of the finite 
system size and can be observed also for static quantities (see [l!4j ). Analogously to statics, finite-size effects can 
be expected to be much less pronounced when looking directly at the wavenumber dependence of a critical quantity. 
Indeed, in Fig. [5}) it is clearly seen that, when approaching the critical point, T smoothly assumes its expected power- 
law oc k z . Far above the critical point, order-parameter modes at low k cross over to the propagating regime and a 
relaxation rate can not be defined anymore. 

In Fig. [St, the effective longitudinal viscosity, vm{k) = (? sR (k) /T(k) = ksT /\T(k)C(k)], computed from the data 
of the relaxation rate and static structure factor, is shown. Far above the critical temperature, vm assumes its 
bare value vi and is practically independent of wavenumber, while at criticality {6 = 0), the expected power-law 
divergence vm oc k~ x is reproduced with reasonable accuracy. Note that the critical enhancement is multiplicative 
and independent of the bare viscosity. Deviations at the lowest k can be at tributed to the relatively strong finite-size 
effects that occur in the static structure factor of the present model |ll4| : at low k, C(k) appears slightly steeper 
than the expected fc _2+I, -power law at the critical point, which is reflected in a wcakcr-than-cxpectcd divergence of 
the longitudinal viscosity. 



2. Shear viscosity 



Turning to the critical behavior of the shear viscosity, we study here only the contribution from the off-diagonal 
component of the thermodynamic pressure tensor, P xy = n(d x p)(d y p), via the Green-Kubo relation, eq. (|68p. Fig. [SJt 
shows the so obtained temperature-dependence of the fluctuation contribution i/ S)Cr it at fc = in the critical region. 
The theoretical prediction of eq. (|68|). represented by the solid curve in Fig. UK, is found to capture the behavior of 
the shear viscosity quite well, after correcting for an overall prefactor of 0(1). To avoid complications associated with 
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FIG. 5: (a) Dependence of the relaxation rate on the reduced temperature 9 = (r c — r)/r c (r is the Landau parameter and r c 
its value at the critical point) for different wavenumbers. The expected power-law decrease of T (dashed line) is rounded off 
when approaching the critical temperature due to the finite system size, (b) Wavenumber-dependence of the relaxation rate 
for different reduced temperatures approaching the critical point from above. Data points in the transition region between 
the overdamped and propagating acoustic regime are omitted. The power-law ~ fc 2- " 7 ~ k corresponds to the van Hove 
prediction for T, where a wavenumber-independent vm is assumed, (c) Effective longitudinal viscosity um (normalized to its 
bare value v{) in dependence of temperature and wavenumber, approaching the critical point from above. Finite-size effects, 
leading to a weaker-than-expected divergence of vm, are noticeable at low k (see text). The lines are drawn as a guide to the 
eye. 



temperature dependent effective critical exponents, mean-field values for 77 and z have been used in the analytical 
evaluation of eq. (|68[) . for simplicity. In Fig. [BL, the wavenumber-dependence of the shear viscosity at the critical 
point obtained from eq. is shown on a double-logarithmic (main plot) and logarithmic-linear (inset) scale. While 
there appears some logarithmic growth at larger k, the plateau at low k suggests that the shear viscosity stays finite 
in 2D, in disagreement with the scaling predictions of eq. (|6"5)) [l30j | . The finiteness of the critical shear viscosity 
can also be seen from its scaling behavior with the system size L (inset to Fig. [SJ;). Interestingly, v s ,cnt is found to 
even decrease with larger L by a power-law with a small exponent of roughly —0.15. From the finite-size behavior of 
the shear-stress correlation function, it is inferred that this behavior arises from both a decrease of the shear-stress 
relaxation time and a decrease of the equal-time autocorrelation of P xy with L (Fig. We find a similar behavior 
also slightly away from the critical point, although the effect is less pronounced there. The decay of the shear-stress 
correlation function is non-exponential and occurs over a shorter timescale than the relaxation of the order parameter. 

A possible reason for the disagreement with the critical scaling predictions of eq. (|68|) might be that, by computing 
the correlation function of the stress tensor P xy in our simulations, we do not properly take into account the renor- 
malization of k, which enters here only as a constant numerical prcfactor. Indeed, if this effect is neglected in eq. (|68l) . 
the scaling exponent changes to z — d — 2i] < in 2D, implying a non-divergent v s , C i-it- An alternative method to 
determine the effective shear viscosity would be to compute the transverse curre nt correlat ion function. However, 
due to the small value of K, which is a constraint of the present LB model (see jll3l Ill4l |). i/ SjC rit remains orders 
of magnitude below its bare value v s . Thus, exceeding computational resources would be necessary to extract the 
fluctuation contribution to the shear viscosity obtained from a correlation function. Further numerical investigation 
of the shear viscosity in a two-dimensional isothermal fluid using alternative simulation methods are thus desirable 
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FIG. 6: (a) Critical fluctuation contribution to the shear viscosity in dependence of the reduced temperature. The data (filled 
symbols) is obtained from simulations using the Green-Kubo relation (|68[) at k = 0. The solid curve represents the theoretical 
prediction of eq. (|68p . (b) Wavenumber-dependence of the critical shear viscosity (at 9 = 0) obtained from eq. (|68p . plotted 
in double- logarithmic representation. For comparsion, the same data is plotted in the inset in log-linear scale, (c) Finite- 
size behavior of the shear stress correlation function (normalized by the system volume V = L 2 ) and the critical fluctuation 
contribution to the shear viscosity, ^ s , C rit (inset). The dashed curve is oc L~ 0,15 . Time is expressed in units of the characteristic 
order-parameter relaxation time 1/T* 

for future work. 



IV. DISCUSSION 

In a conventional single-component fluid (model H) , the dominant transport mechanism is heat diffusion, while sound 
waves are decoupled from the order-parameter dynamics [H 1(1 l54l. f&sj] . In contrast, under isothermal conditions, heat 
diffusion is absent and order-parameter fluctuations can relax only via sound waves. Using a mode-coupling approach, 
it has been argued here that, below four dimensions, the order-parameter dynamics of the critical isothermal fluid 
is characterized by model-A-type behavior. This means that, at long wavelengths, sound waves are ovcrdamped due 
to a strongly diverging bulk viscosity, V\ oc £ x . The relaxation rate, T = c 2 s /vi, scales as with a dynamic critical 
exponent of z = j/i> + x = 2 — r/ + x, where j/v represents the contribution from the isothermal speed of sound 
and x ~ 1.77? represents the "non-classical" contribution from the renormalized bulk viscosity. While in model H, 
the divergence of the kinetic coefficient arises due to the advective coupling to the transverse velocity modes @, H3] , 
it is due to the nonlinear thermodynamic pressure in the isothermal case. Longitudinal and transverse currents are 
approximately decoupled for d close to 4 - a property that becomes exact in the linear case. In 2D, our simulations of 
the fluctuating hydrodynamic equations yield a value ofz w 2.2 ±0.1 and x w 0.45 ±0.1, in reasonable agreement with 
theoretical predictions and Monte-Carlo simulations of model A. The scaling theory predicts the shear viscosity to 
remain finite for d > 2 and weakly diverge by a power-law in two dimensions. This divergence could, however, not be 
observed within the present simulation approach. The essential differences in the critical dynamics of an isothermal 
and an ordinary fluid (model H) are collected in Tab. |U 

It is interesting to compare the present findings also to the situation in hydrodynamic models of the glass transition 
[35I [111 [39l |95| |. where the isothermal compressible Navier-Stokes equations have been investigated in conjunction 
with a purely Gaussian free energy. There, the density correlation function shows an an anomalously slow decay at 
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low temperatures, accompanied by a strong increase of the bulk viscosity. Quite analogously to the case in critical 
dynamics, this is a generic mode-coupling effect caused by a nonlinear pressure term. In the case of a supercooled 
liquid, however, the dominant contribution arises from the quadratic pressure nonlinearity [r</> 2 in cq. (|20|)]. whereas 
in the critical fluid, this term turns out to be irrelevant due to the smallness of r. 

The present model is particularly interesting in the two-dimensional case, where a possible experimental realization 
of the isothermal condition might be achievable. While the present scaling considerations indicate the divergence 
of additional contributions to the bulk viscosity beyond the model-A term, our simulation results suggest that the 
model- A- type critical behavior essentially persists also in 2D, with possible corrections to critical exponents being 
small, at least. For future work, it will be interesting to treat the isothermal fluid model within a rcnormalization- 
group approach and also derive more detailed predictions in the two dimensional case. In order to clarify the critical 
behavior of the shear viscosity in 2D, alternative simulation methods could be used. 

The predictions obtained in this work might be experimentally testable on single-component monolayer films that 
admit for liquid- vapor-like phase-separation below a critical point (Til [l2j. Of course, t he presen t model is highly 
idealized in that it neglects the possible influence of electrost atic l o ng-r ange interactions jl20l - ll23| . friction between 
fluid and substrate, and hydrodynamic back-coupling [H, [88|, Il24l . Il25| . Also, it is assumed that the rate of heat 
transfer between fluid and s ubstr ate is sufficiently large to provide an effective isothermal environment for the critical 
fluctuations (see, however, |l26| |). Since the long- wavelength dynamics of a fluid becomes arbitrarily slow upon 
approaching the critical point, one might expect that even a relatively small thermal coupling will actually be sufficient. 

In general, crossover behavior between different dynamic universality classes of a single-component fluid film is 
expected: in case of negligible friction, model H is obtained for vanishing thermal coupling and, as shown in this 
work, model A for perfect thermal coupling. In the opp osite case of large friction, it is expected that model B (i.e., a 
purely diffusive order-parameter transport) results j79l [88j. 
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Appendix A: Further contributions to the self-energy 



The contributions to the order-parameter self-energy at two-loop order are diagrammatically represented as: 




Here, the unlabeled dots represent the possible couplings r, K, vi, vt associated with the three-point vertices. Open 
circles and arrows (representing correlation or response functions) are understood to be present on some of the 
internal lines. Also, some diagrams exist with the solid lines replaced by wavy lines representing the transverse 
current response/correlation functions. Regarding the scaling behavior, however, these modified diagrams need not 
be explicitly evaluated, since all three-point vertices scale in the same way. Taking into account the proper critical 
behavior of the couplings, the dynamic part of each diagram is found to scale oc £ 2_d (up to corrections of exponents 
of 0(j])) and thus gives negligible contributions to the renormalized viscosity for d > 2. 
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relation like eq. (|3ip . Due to the presence of the 1/p-nonlinearity, however, there is no simple FDT like eq. (|30[1 connecting 
the full correlation and response function in the compressible fluid, except in the small-fc limit [36], (3S| . An FDT can be 
proven for the nonlinear oscillator eq. (J26J if the l/p-nonlinearity is neglected, see [9j| 
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correlation function decays non-exponentially and the dynamics deviates from pure model-A behavior. 

[130] Although the shear- viscosity at k — when plotted against reduced temperature (Fig. [6^.) exhibits an extended plateau 
towards 8 — y as well, this effect can not be unambiguously attributed to the non-divergent nature of f s , cr it, as finite-size 
effects are expected to contribute significantly here (cf. Fig. [5^,b, where a similar effect is seen for the relaxation rate). 



